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Abstract 

The experimental data on p, A, f2“-baryons and the corresponding antibaryons 
spectra obtained by different collaborations are compared with the results of the calcu¬ 
lations performed into the frame of the Quark-Gluon String Model. The contribution 
of String Junction diffusion and the inelastic screening corrections are accounted for 
in the theoretical calculations. The predictions of the Quark-Gluon String Model both 
for pp and pA collisions are extended up to the LHG energies. 
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1 Introduction 


The study of hadron production has a long history in high-energy particle and nuclear 
physics, as well as in cosmic-ray physics. The absolute yields and the transverse mo¬ 
mentum (pt) spectra of identihed hadrons in high-energy hadron-hadron collisions are 
among the most basic physical observables. They can be used to test the predictions 
for non-perturbative quantum chromodynamics (QCD) processes like hadronization 
and soft-parton interactions, and the validity of their corresponding implementation in 
Monte Carlo (MC) event generators [U [21 E] • 

Among the different available phenomenological models, in this paper we will con¬ 
sider the Quark-Gluon String Model (QGSM) [H E], based on the Dual Topological 
Unitarization, Regge phenomenology, and nonperturbative notions of QCD, and that 
has been successfully used for the description of multiparticle production processes in 
hadron-hadron [HlElElISlllo] and hadron-nucleus HH El na [H E] collisions. 

In the QGSM high energy interactions are considered as proceeding via the exchange 
of one or several Pomerons, and all elastic and inelastic processes result from cutting 
through or between Pomerons [16]. Inclusive spectra of hadrons are related to the 
corresponding fragmentation functions of quarks and diquarks, which are constructed 
using the Reggeon counting rules na. 

In the case of interaction with a nuclear target, the Multiple Scattering Theory 
(Gribov-Glauber Theory) is used, what allows to consider the interaction with nucleus 
as the superposition of interactions with different numbers of target nucleons m- 

The QGSM gives a reasonable description of a lot of experimental data on baryons 
and strange hyperons production in hadron-hadron collisions for wide energy range 
going from hxed target experiments up to LHG [SI [ZllSl US] [9| HH |l2l [131 HH [15] . In 
this paper we extend the model to the study the multistrange hyperon production. 

The multistrange baryons, E~ (dss), and (sss), are particulary important in 
high energy particle and nuclear physics, due to their dominant strange quark (s- 
quark) content. Since the initial-state colliding projectiles contain no strange valence 
quarks, all particles in the hnal state with non-zero strangeness quantum number should 
have been created in the course of the collision. Moreover, the energy of the Large 
Hadron Gollider (LHG) and its high luminosity allow for an abundant production of 
strange hadrons. These two factors make multi-strange baryons a valuable probe in 
understanding the particle production mechanisms in high energy collisions. 

A remarkable feature of strangeness production is that the production of each ad¬ 
ditional strange quark featuring in the secondary baryons, i.e., the production rate of 
secondary B{qqs) over secondary B{qqq), then of B{qss) over B{qqs), and, hnally, of 
B{qss) over B{sss), is affected by one universal strangeness suppression factor. A*: 

^ Bjqqs) ^ B{qss) ^ B{qss) 

B{qqq) B{qqs) B{sss) 
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together with some simple combinatorics 

It will be shown below that A* slightly increases with initial energy, the experimental 
data favoring the value As=0.22 at hxed target eneries < 30 GeV), that becomes 

As=0.32 at larger RHIC and LHC energies. This energy dependence can be connected 
to an increase of hardness of the interacton, i.e. to the increase of the secondaries 
average transverse momenta IZD]. 

Signihcant differences in the yields of baryons and antibaryons in the central (midra¬ 
pidity) region are evident at not very high energies. This effect can be naturally ex¬ 
plained P ini Uni EH EH E3] by the structure of baryons in QGSM, i.e. baryons of three 
valence quarks together with a specihc conhguration of the gluon held, called String 
Junction [211 [251 l&l [271. 

At very high energies the contribution of the enhanced Reggeon diagrams becomes 
important, leading to a new phenomenological effect, the suppression of the inclusive 
density of secondaries [2H] in the central (midrapidity) region. 

We present a brief description of the model in Section 2. In Section 3 we show 
the comparison of our numerical calculations with the existing experimental data on 
secondary baryons and antibaryons production in pp and pA collisions. Gonclusion 
remarks are given in Section 4. 


2 QGSM formalism 

2.1 General approach 

The QGSM [H El [7] allows one to make quantitative predictions for different features 
of multiparticle production, in particular, for the inclusive densities of different secon¬ 
daries, both in the central and in the beam fragmentation regions. 

In QGSM, each exchanged Pomeron corresponds to a cylindrical diagram, and thus, 
when cutting one Pomeron, two showers of secondaries are produced (see Fig. 1 a,b). 

The inclusive spectrum of a secondary hadron h is then determined by the convolu¬ 
tion of the diquark, valence quark, and sea quark distributions, u{x, n), in the incident 
particles, with the fragmentation functions, G^{z), of quarks and diquarks into the 
secondary hadron h. Here n is the number of cutted Pomerons in the considered dia¬ 
grams. Both the distributions and the fragmentation functions are constructed using 
the Reggeon counting rules nzi. 

In particular, in the case of n > 1, i.e. in the case of multipomeron exchange, the 
distributions of valence quarks and diquarks are softened due to the appearance of a 
sea quark contribution. 

The details of the model are presented in HlEllTl [in]. The average number of 
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Figure 1: (a) Cylindrical diagram representing a Pomeron exchange within the DTU classification 
(quarks are shown by solid lines); (b) One cut of the cylindrical diagram corresponding to the single- 
Pomeron exchange contribution in inelastic pp scattering. 


exchanged Pomerons {n)pp slowly increases with the energy. The values of the Pomeron 
parameters have been taken from [7]. 

For a nucleon target, the inclusive rapidity, y, or Feynman-x, xp, spectrum of a 
secondary hadron h has the form |1]: 


dn 

dy 


xe da 
^inel dx E 


oo 

= Y,Wn- , 


n=l 


( 2 ) 


where the functions (p^ix) determine the contribution of diagrams with n cut Pomerons, 
Wn ID is the relative weight of this diagram. Here we neglect the numerically small 
contribution of diffraction dissociation processes. 

In the case of pp collisions: 

0nE) = fqq{x+,n)-fl^{x-,n)+f’P^{x+,n)-fl^q{x-,n)+2{n-l)-f’^{x+,n)-f’^{x-,n) , (3) 


x± 


1 

2 


[\jAm\/s + ± x] , 


(4) 


where fqq, fq, and fs correspond to the contributions of diquarks, valence quarks, and 
sea quarks, respectively. 


These contributions are determined by the convolution of the diquark and quark 
distributions with the fragmentation functions, e.g.. 


fg{x+,n) = f Uq{xi,n) ■ G'^{x+/xi)dxi 

J 


(5) 
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2.2 String Junction contribution 


In the string models, baryons are considered as configurations consisting of three con¬ 
nected strings (related to three valence quarks), called String Junction (SJ) 125112^ 
I2Z], this picture leading to some quite general phenomenological predictions. 

The production of a baryon-antibaryon pair in the central region usually occurs via 
SJ-SJ pair production (SJ has upper color indices, whereas anti-SJ {SJ) has lower 
indices), which then combines with sea quarks and sea antiquarks into a BB pair [26], 
as it is shown in Fig. 2a. 

d g 
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Figure 2: QGSM diagrams describing secondary baryon production: (a) usual BB central production 
with production of new SF pair; (b) initial SF together with two valence quarks and one sea quark; 
(c) initial SF together with one valence quark and two sea quarks; (d) initial SF together with three 
sea quarks. 

In the case of central production of baryon—antibaryon pairs (Fig. 2a) the fragmen¬ 
tation functions of any quarks and diquarks into antibaryons are proportional to 

( 6 ) 

where parameters and take different values for different antibaryon states. The 
parameter (3^ was described in [101 HB 1131 EO]. In the central region, z —)■ 0, the 
values of the parameter in the fragmentation functions of Eq. 5 for p. A, S, and 
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productions are connected by the following relation |13j : 


Oj, : OJ : o, : Ojj = 1 : ^(5/2) • A. : v'(3/4) ■ : v'(l/4) ■ Af , (7) 

where is the strangeness suppression factor (see Introduction), and the numerical 
factors ^5/2, ^3/4, and y^l/4 are determined by quark combinatorics (THl HSl EH]- 
The value of the parameter ap for p production was determined by comparison with 
the experimental data on antiproton production in pp collisions at different energies to 

be dH ESI [29] 

Op = 0.18. (8) 

The comparatively large value for A production comes from the fact that many 
strange hyperon resonances (e.g. S*(1385)) decay into A and not into S. 

By using the relations between ap, a-^, and a-^ in Eq. ([^, that derive from 
simple quark combinatorics [101 EH ES], we can calculate the constants in Eq. ([^ 
corresponding to the central production of any strange and multistrange antibaryon in 
the diagrams describing the central BB pair production (see Fig. 2a). 

However, in the processes with incident baryons another possibility exists to produce 
a secondary baryon in the central region, called SJ diffusion (see hgs. 2b,c,d). 

To obtain the net baryon charge, and according to ref. m, we consider three 
different possibilities. The hrst one is the fragmentation of the diquark giving rise to 
a leading baryon (Fig. 2b). A second possibility is to produce a leading meson in the 
hrst break-up of the string and a baryon in a subsequent break-up |T^ (Fig. 2c). In 
these two hrst cases the baryon number transfer is possible only along short distances 
in rapidity. In the third case, shown in Fig. 2d, both initial valence quarks recombine 
with sea antiquarks into mesons, M, while a secondary baryon is formed by the SJ 
together with three sea quarks PEo]. 

The fragmentation functions for the secondary baryon B production corresponding 
to the three processes shown in Fig. 2b, 2c, and 2d can be written as follows (see [101 
for more details): 







(9) 

Gfsi 


- ttp 

■vfs ■ 


(10) 

Gfsi 



■ e ■ V 

P . . (1 _ ^ 

(11) 


where 

Up = 1.33 (12) 

is the normalization parameter, obtained earlier from the description of proton spectra 
in pp collisions [niE 3 l| 29 ], and Vqg, Vgg, are the relative probabilities for diherent 
baryons production that can be found by simple quark combinatorics [THl EH] (see 
Table 1), and the factor e accounts for the small probability of the process shown in 
Fig. 2d to occur with respect to the processes in Figs. 2b,c. 
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Table 1: The quark combinatorics providing the values of the parameter vf in eqs. (8—10). Here 
S/L — \g and 2L S = 1. 

The probabilities depend on the vaule of the strangeness suppression 

factor As- In [TnilI3] we used As=0.32. In the present paper we show that A^ = 0.22 is 
in a better agreement with the main sample of the data at intermediate energies (see 
section 3). 

The fraction 2 ; of the incident baryon energy carried by the secondary baryon de¬ 
creases from Fig. 2b to Fig. 2d. Only the processes in Fig. 2d can contribute to the 
inclusive spectra in the central region at high energies if the value of the intercept of 
the SJ exchanged Regge-trajectory, asj, is large enough. The analysis in [23] gives a 
value of asj = 0.5 ± 0.1, that is in agreement with the ALICE Collaboration result, 

~ 0.5 [30]. In the present calculations we use [23] : 

a^j = 0.5 and £ = 0.0757. (13) 

2.3 Interaction with nuclei at high energies and inelastic screen¬ 
ing effects 

In the calculation of the inclusive spectra of secondaries produced in pA collisions we 
should consider the possibility of one or several Pomeron cuts in each of the f blobs 
of proton-nucleon inelastic interactions. For example, in Fig. 3 it is shown one of the 
diagrams contributing to the inelastic interaction of a beam proton with two target 
nucleons. In the blob of the proton-nucleonl interaction one Pomeron is cut, while in 
the blob of the proton-nucleon2 interaction two Pomerons are cut. 

The contribution of the diagram in Fig. 3 to the inclusive spectrum is: 

2 ■ ^v(2) ■ ■ {/‘,q+, 3) ■ 1) + 

3) • fgg{x-, 1) -b f^{x+, 3) • [fqq{x-, 2) -f 2) + 

2-f^{x.,2)]} , (14) 
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Figure 3: One of the diagrams for the inelastic interaction of one incident proton with two target 
nucleons and N 2 in a pA collision. 

where WpA{2) is the probability of interaction with namely two target nucleons. 

It is essential to take into account all diagrams with every possible Pomeron con- 
hguration and its corresponding permutations. The diquark and quark distributions 
and the fragmentation functions are the same as in the case of pN interaction. 

The total number of exchanged Pomerons becomes as large as 

{n)pA ~ {j^)pA ■ {n)pN , (15) 

where {i')pA is the average number of inelastic collisions inside the nucleus (about 4 for 
heavy nuclei at hxed target energies). 

The process shown in Fig. 3 satishes [SH |32l |33l |3l] the condition that the absorp¬ 
tive parts of the hadron-nucleus amplitude are determined by the combination of the 
absorptive parts of the hadron-nucleon amplitudes. 

The QGSM assumption of the superposition picture of hadron-nucleus interactions 
gives a reasonable description ilElES] of the inclusive spectra of different secondaries 
produced at energies y/sNN = 14—30 GeV. 

At RHIG energies the situation drastically changes. The spectra of secondaries pro¬ 
duced in pp collisions are described rather well [9] , but the RHIG experimental data for 
Au + Au collisions [361137] give clear evidence of the inclusive density suppression effects 
which reduce by a factor ~0.5 the midrapidity inclusive density, when compared to the 
predictions based on the superposition picture |2Hl EH] • This reduction can be explained 
by the inelastic screening corrections connected to multipomeron interactions [2H| (see 
Fig. 4). 

At energies ^/s]^ < 30—40 GeV, the inelastic processes are determined by the 
production of one (Fig. 4a) or several (Fig. 4c) multiperipheral ladders, and the corre- 
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Figure 4: (a) Multiperipheral ladder corresponding to the inclusive cross section of diagram (b), and 
(c) fusion of several ladders corresponding to the inclusive cross section of diagram (d). 


spending inclusive cross sections are described by the diagrams of figs. 4b and 4d. 

In accordance with the Parton Model nnun], the fusion of multiperipheral ladders 
shown in Fig. 4c becomes more and more important with the increase of the energy, 
resulting in the reduction of the inclusive density of secondaries. Such processes corre¬ 
spond to the enhancement Reggeon diagrams of the type of Fig. 4d, and to even more 
complicate ones. All these diagrams are proportional to the squared longitudinal form 
factor of nucleus A |2H]. On the contrary, the contribution of such diagrams becomes 
negligible when the energy decreases, and they are also negligible for pp collisions up 
to LHC energies. 

In any case, all quantitative estimations of the importance of the contribution of the 
enhanced diagrams are model dependent, since the numerical weight of the contribution 
of the multipomeron diagrams is rather unclear due to the many unknown vertices in 
these diagrams. In some models the number of unknown parameters can be reduced, 
and, as an example, in reference [28] the Schwimmer model ns was used for the 
numerical estimations. 

Another possibility to weight the contribution of the diagrams with Pomeron inter¬ 
action comes gal m ig SSI ST] from Percolation Theory. The percolation approach 
and its previous version, the String Fusion Model gUSnilSQl, predicted the multiplicity 
suppression seen at RHIC energies, long before any RHIC data were measured. 

In order to account for the percolation (inelastic screening) effects in the QGSM, 
it is technically more simple [22] to consider a maximal number of Pomerons n^ax 
emitted by one nucleon in the central region that can be cut. These cut Pomerons lead 
to the different hnal states. Then, the contributions of all diagrams with n < Umax are 
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accounted for as at lower energies. The larger number of Pomerons n > Umax can also 
be emitted obeying the unitarity constraint, but due to the fusion in the hnal state (at 
the quark-gluon string stage), the cut of n > Umax Pomerons results in the same final 
state as the cut of Umax Pomerons. 

By doing this, all model calculations become very similar to those in the percolation 
approach. The QGSM fragmentation formalism allows one to calculate the integrated 
over pt spectra of different secondaries as functions of the rapidity and xp- 

In this frame, we obtain a reasonable agreement with the experimental data on 
the inclusive spectra of secondaries produced in d -|- Au collisions at RHIC energy [22] 
with a value Umax = 13, and in p -|- Pb collisions at LHC energy |13] with the value 
^max = 23. It has been shown in [5T] that the number of strings that can be used for 
the secondary production should increase with the initial energy. 


3 Numerical results 

3.1 The pp collisions 

The existing experimental data for the integrated over the whole range of pp energy 
dependence of p, p. A, A, and hyperons production density dn/dy{\ y \< 0.5) 

in pp collisions, from hxed target to LHC energies [521 ESI EH EH EH EH ESI EH EH Elj 
are compared with the results of the QCSM calculations in Table 2 and in Fig. 5. 

The QCSM calculation for secondary p and p production in pp collisions performed 
with the same values of the parameters ap=0.18 (see Eq. 7) and 0^=1.33 (see Eq. 11) 
used in our previous papers [mEH, are in reasonable agreement with the experimental 
data [SHEHEHEIEHEHEZ]. 

For all baryons (p. A, S, C) and for their corresponding antibaryons, the calculated 
values for dn/dy(| y \< 0.5) are shown in Fig. 5 by full curves, and for antibaryons by 
dashed curves. The theoretical curves and the data on p and p production in Fig. 5 
are shown as multiplied by a factor 10. 

The data on A and A [53l [52] productions are reasonably described at the lower en¬ 
ergies with the value Xs=0.22, that should be increased up to As=0.32 at ^/s >60 CeV. 
This increase of the value of can be connected to the increase in the hardness of 
the interaction, since then the non-zero value of the s-quark mass becomes not so sig- 
nihcant. Possibly, the same effect also leads to the increase of average pp [2T]. At 
intermediate energies the hardness of the interaction is relatively small and the pro¬ 
duction of the s-quark, with a mass rris ~ 0.3GeIA, is seriously suppressed. At higher 
enetgies, the hardness of the interaction increases and the effect of the hnite s-quark 
mass becomes not so important, what is refleted by the effective increase of the value 
of the As parameter. At extremely high energies, as the hardness of the interaction is 


10 



Figure 5: The energy dependence of baryon and antibaryon production inclusive densities 
dn/dyd y \< 0.5) in the midrapidity region in pp collisions (baryons are shown by closed points 
and full curves, and antibaryons by open points and dashed curves). 


much larger, we expect —?• 1. 

The yields of and produced at 158 GeV/c [M] and at LHC energies of 0.9 
and 7 Tev [53] are well reproduced by correspondingly using the same values of Xs=0.22 
and As=0.32. 

In Fig. 5 we also present our prediction for and hyperon production. 

We do not include in our analysis the ALICE Collaboration data for 
and C hyperons production [69] at 7 TeV for the reason that these data were measured 
at transverse momenta pp > 0.6 GeV/c, and therefore they can not be compared to 
integrated over the whole range of pp rapidity distributions. 
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ys (GeV) 

Reaction 

QGSM 

Experiment dn/dy (| y \< 0.5) 

17.29 

P + P -tp 

0.0605 

0.0736 ±0.015 [M] 


P + P -tp 

0.024 

0.054 ±0.0014 [Ml 

27.45 

P + P 

0.0576 

0.0736 ±0.015 |M] 


P + P 

0.032 

0.01869 ±0.001 [M] 

62.4 

P + P 

0.064 

0.040 ± 0.003 [63! 


P + P 

0.0477 

0.022 ±0.002 [63! 

200. 

P + P -tp 

0.0454 

0.044 ± 0.004 [Ml 


P + P -tp 

0.0392 

0.037 ±0.003 [Ml 

900. 

P + P ^P 

0.111 

0.106 ±0.005 IMj 


P + P ^P 

0.111 

0.083 ± 0.006 [Ml 


P + P ^P 

0.107 

0.079 ± 0.005 [Ml 


P + P ^P 

0.107 

0.104 ±0.003 |65| 

2760. 

P + P -tp 

0.138 

0.124 ±0.003 |M] 


P + P -tp 

0.135 

0.1235 ±0.0035 [M] 

7000. 

P + P -tp 

0.161 

0.15 ±0.003 |65| 



0.161 

0.124 ±0.009 IMj 


P + P ^P 

0.160 

0.147 ± 0.0035 [M] 



0.160 

0.123 ±0.01 [M] 

17.29 

p + p A 

0.0219 

0.017 ±0.0004 IM] 


p + p A 

0.0065 

0.0057 ±0.0002 [M] 

23.79 

p + p A 

0.0227 

0.0132 ± 0.003 [57] 


p + p A 

0.0105 

0.010 ±0.003 [Ml 

27.45 

p + p ^ A 

0.0233 

0.0091 ±0.0015 [60] 


p + p ^ A 

0.0116 

0.017 ±0.003 [M] 

200. 

p + p ^ A 

0.0454 

0.0436 ± 0.0008 ± 0.004 [M] 


p + p A 

0.0392 

0.0398 ± 0.0008 ± 0.0037 [M] 

900. 

p + p A 

0.0659 

0.26 ±0.01 [52] 

7000. 

p + p A 

0.0987 

0.27 ±0.01 [M] 

900. 

p + p ^ A 

0.0659 

0.108 ±0.001 ±0.012 [531 

7000. 

p + p A 

0.0987 

0.189 ±0.001 ± 0.022 [M] 

17.29 

p + p ^ s- 

0.000839 

0.000799 ± 0.000095 [54] 


p + p ^ 

0.000236 

0.000381 ± 0.000057 [M] 

200. 

p + p ^ s- 

0.00542 

0.0026 ± 0.0002 ± 0.0009 [M] 


p + p ^ S 

0.00487 

0.0029 ± 0.0003 ± 0.006 [M] 

900. 

p + p ^ s- 

0.00942 

0.011 ±0.001 ±0.001 [M] 

7000. 

P + P ^ S” 

0.0159 

0.021 ±0.001 ±0.003 [53| 


P + P ^ S- 

0.0159 

0.008 ±0.001 ±0.007 [69] 


p + p s’*" 

0.0158 

0.0078 ± 0.001 ± 0.007 [M] 

200. 

p -l- p — y ri -|- ^2”^ 

0.0000634 

0.00034 ± 0.00016 ± 0.0005 [M] 

7000. 

p + p — >■ 

0.00132 

0.00067 ± 0.00003 ± 0.00008 [M] 


p + p — >■ 

0.00134 

0.00068 ± 0.00003 ± 0.00008 [69] 


Table 2. The experimental data on the inclusive density dn/dy of baryons and antibaryons produced 
in pp collisions at different energies and their corresponding description by the QGSM. 

In Fig. 6 we compare our QGSM calculations with the experimental data on the 
rapidity dependence of the production cross section da/dy, integrated over the whole pp 
range, of proton and antiproton production in pp collisions at two energies, 158 GeV/c [M] 
and 400 Gev/c [62]. The full curve corresponds to the proton production and the 
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dashed curve to p production, both at 158 GeV/c [5l], while the dashed-dotted curve 
corresponds to p production and the dotted curve to p production at 400 GeV/c [62] . 

On the lower panel of Fig. 6, the A and A production cross section da/dy, at 
158 GeV/c [51| and at 405 GeV/c (HUj in pp collisions are compared with the results of 
the QGSM calculations. The full curve corresponds to A production at 158 GeV/c and 
the dashed curve to A production at the same energy, while the dashed-dotted curve 
corresponds to A production and the dotted curve to A production, both at 405 Gev/c. 
All curves for A and A production have been obtained with a value of As=0.22. As one 
can see in Fig. 6, the model calculations agree with the experimetal data. 
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Figure 6: The QGSM results for (upper panel) the rapidity dependence of the production densi¬ 
ties dn/dy of p and p at 158 GeV/c and 400 GeV/c, and their comparison with the experimental 
data |54ll62|, and for (lower panel) the A and A cross sections at 158 GeV/c and 405 GeV/c, and their 
comparison with the experimental data |54[ I60j (see the main text por the description of the different 
curves). 
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In Fig. 7 we compare the QGSM calculations for the rapidity dependence of the 
inclusive cross section da jdy integrated over the whole pt range of and produced 
in pp collisions at 158 GeV/c [M]- Also here, the theoretical curves have been data 
obtained by using the vane As=0.22. The full curve corresponds to production and 
the dashed line to the production of 

On the upper panel of Fig. 8 we show the data on the inclusive density {l/nev)dn/dy 
for A [521 El] and [53] hyperons production at both 0.9 and 7 TeV. Here it has to 
be noticed that the data by ATLAS [52] and GMS [53] collaborations at LHG energies 
0.9 and 7 TeV are inconsistent with each other, and, in particular, the ATLAS data 
show a weak energy dependence, while the GMS data would indicate a significantly 
stronger energy dependence. 

As one can see in Fig. 8, the QGSM calculations are in a satisfactory agreement for 
the cases of A production at 0.9 TeV and for production, both at 0.9 and 7 TeV, 
and they are below the experimental points for A production at 7 TeV. 



y 

Figure 7: The QGSM results for the rapidity dependence of the inclusive cross section da/dy of 
and s'*" productions in pp collisions at 158 GeV/c, and their comparison with the experimental 
data [Mj. The full curve corresponds to S“ and the dashed curve to S^ production. 

On the lower panel of Fig. 8 we show the experimental data on rapidity dependence 
of ratios of A/S“ densities at 0.9 and 7 TeV [53] and their comparison with the corre¬ 
sponding QGSM results. The QGSM reasonably reproduces both the values and the 
rapidity dependences of the A/S“ ratios at the two energies. 
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Figure 8: The rapidity dependence of (upper panel) the inclusive density {\/nf,v)dn/dy of A and 
hyperons production in the midrapidity region for pp collisions at LHC energies |53l 152) , and of (lower 
panel) of the ratio of yields of A/S“ hyperons at LHC energies |53j in the midrapidity region for pp 
collisions [5S]. The corresponding QGSM predictions are shown by the full curve (-^/s = 7 Tev) and 
by the dashed curve [^/s = 0.9 TeV). 


3.2 The pA collisions 

To perform the calculations corresponding to the case of proton-nucleus interactions we 
have used the same values of the parameters Op, a^, and as for pp collisions, and we 
also keep nmax= 23 (see subsection 2.3) to account for the inelastic nuclear screening 
effects. 

In Fig. 9 and Table 3 we compare the results of the QGSM calculations for the 
A-dependences of p, p, and A, A, S'*", and G"'' hyperons produced on different 
nucleus targets with the corresponding experimental data IZH Ca [731 [H ES]. The 
QGSM results at 158 GeV/c are presented for all baryons. For A and A production 
two sets of measurements at 158 Gev/c and at 200 GeV/c exist. At 158 GeV/c, the 
full curves correspond to the production of baryons and the dashed curves to that of 
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Figure 9: Comparison of experimental data on the A-dependence of midrapidity density dn/dy 
(I y |< 0.5) of p, p, and A, A, S , and ft hyperons produced at 158 Gev/c and 200 
GeV/c [Til |73l [ZH [75] , with the results of corresponding QGSM calculations (see the main text por 
the description of the different curves). 


antibaryons. The curves corresponding to the calculations at 200 GeV/c are very close 
to those obtained at 158 GeV/c, and then not shown. The QGSM curves are here in 
a reasonable agreement with the experimental data. 

The experimental data on the rapidity density dn/dy{\ y |< 0.5) for A, A, 

G", and hyperons productions at 158 GeV/c and for p and p production at 
^/s=5.02 TeV [70] in pPb collisions |7I1 [73] [TH [75] are presented in Fig. 10, together 
with the corresponding QGSM predictions from hxed targets up to LHG energies. 

Here, as in Fig. 5, the curves for p and p production and the GMS Gollaboration ex¬ 
perimental data on average [p + p)/2 production inp+Ph collisions at ^/s=5.02 TeV [TO] , 
have been multiplied by a factor 10 to make them distinct from the A and A case. 

In Fig. 11 we compare the experimental data on the rapidity dependences of the 
denstity dn/dy of p and p produced in pC collisions measured by the NA49 Gollabora¬ 
tion at 158 GeV/c [72], together with the corresponding QGSM calculations. The full 
curve corresponds to p production and the dashed one to p production. The agreement 
between the QGSM results and the experimental data shown in the figure is rather 
good. 
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vG (GeV) 

Reaction 

QGSM 

Experiment dn/dy(| y \< 0.5) 

17.2 

p -1- G p 

0.0720 

0.12157 ±0.005 


p -1- G p 

0.0285 

0.02656 ±0.0011 [72] 

5020. 

p -1- Pb (p -|-p)/2 

0.440 

0.5 ±0.05 |7D| 

17.2 

p -t Be — )■ A 

0.0241 

0.0334 ± 0.0005 ± 0.003 |7T] 


p -t Be — )■ A 

0.0107 

0.011 ±0.0002 ±0.001 ini 

19.42 (200 GeV/c) 

p -t Ar A 

0.0320 

0.02 ±0.01 [73] 

19.42 (200 GeV/c) 

P -t Xe ^ A 

0.0413 

0.03 ±0.015 [zg 

19.42 (200 GeV/c) 

p -1- Au A 

0.0450 

0.05 ±0.01 izg 

19.42 (200 GeV/c) 

p -1- Au A 

0.0193 

0.0095 ±0.01 [III 

19.42 (200 GeV/c) 

p -t S — >■ A 

0.0306 

0.06 ± 0.008 [75| 

19.42 (200 GeV/c) 

p -t S — >■ A 

0.0131 

0.015 ±0.003 US] 

17.2 

p -t- Pb A 

0.0458 

0.060 ± 0.002 ±0.006 [7T| 

(m.b.) 

p -t Pb A 

0.0173 

0.015 ±0.001 ±0.002 [n] 

17.2 

p -1- Be — 

0.0011 

0.0015 ± 0.0001 ± 0.0002 HI] 

(m.b.) 

p -1- Be — s'*" 

0.000482 

0.0007 ± 0.0001 ± 0.0002 |7l] 

17.2 

p -1- Pb ^ S" 

0.00217 

0.0030 ± 0.0002 ± 0.0003 [Hj 

(m.b.) 

p -b Pb ^ 

0.000726 

0.0012 ± 0.0001 ± 0.0001 HI] 

17.2 

p -l- Be — y 

0.0000434 

0.00012 ± 0.00006 ± 0.00002 HI] 

(m.b.) 

p -b Be —>■ 

0.0000155 

0.00004 ± 0.00002 ± 0.00001 HI] 

17.2 

p -b Pb —>■ 

0.0000749 

0.00022 ± 0.00008 ± 0.00003 HI] 

(m.b.) 

p -b Pb —>■ 

0.0000195 

0.00005 ± 0.00003 ± 0.00002 HI] 


Table 3. Experimental data for baryons and antibaryons production densities dn/dy (| ?/ |< 0.5) in 
proton-nucleus collisions at different energies, and the corresponding description by the QGSM. 

In Fig. 12, the rapidity dependence of the production denstity of A hyperon in 
proton collisions on Ar and Xe nuclei [73] (upper panel) and of A and A hyperons 
produced in proton collisions on S and Au nuclei [HI [75l [76| (lower panel), all at 
200 GeV/c, are compared with the results of the corresponding QGSM calculations. 
The full curve on the upper panel corresponds to A production on Ar nucleus and the 
dash-dotted curve to A production on Xe nucleus. On the lower panel, the full curve 
corresponds to A production on Au nucleus, and the dash-dotted curve to A production 
on S nucleus. Also on the lower panel, the dashed curve corresponds to A production 
on Au nucleus and the dotted curve to A production on S nucleus. 

One can generally affirm that QGSM results on the rapidity dependence of the 
production densities dn/dy of hyperons in pA collisions are in a good agreement with 
the experimental data at 200 GeV/c, as it is shown in hgs. 11 and 12. 
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Figure 10: The energy dependence of the baryon production density dn/dy(\ y \< 0.5) in the 
midrapidity region for pPb collisions. The experimental data for average proton production {p + pl2) 
at 5.02 TeV are by the CMS Collaboration ITD]. The experimental data on hyperon production were 
measured by the NA57 Collaboration at 158 GeV/c [7T]. The curves are the same as in Fig. 4. 

4 Conclusion 

The QGSM, based on the Dual Topological Unitarization, Regge phenomenology, and 
nonperturbative notions of QCD, provides a reasonable description of strange and 
multistrange hyperons production, as well as of the production of their corresponding 
antiparticles, both in pp and pA collisions for a wide range of energies, by including into 
the analysis the contribution of String Junction diffusion and the inelastic screening 
corrections. 

It is important to note that in the QGSM the extension of the calculations of the 
hyperon (antihyperon) productions from pp collisions to proton-nucleus collisions does 
not require any additional model parameters. The accuracy of our calculations can be 
estimated to be on the level of ~ 15%. 

The values of the parameters Op and in the model obtained from previous calcula¬ 
tions at lower energies are now conhrmed by the good description of p and p production 
at higher energies. The estimation of the production of all hyperons and antihyperons 
is governed by the value of only one parameter A^, which is the same for every unit 
of strangeness, but that it depends on energy from the value As=0.22 at hxed target 
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Figure 11: The rapidity dependence of the production denstities of p and p baryons produced in 
pC collisions measured by the NA49 Collaboration at 158 GeV/c [72], together with the results of the 
QGSM calculations. The theoretical curves correspond to those in Fig 4. 

energies, up to As=0.32 at LHC energies. In the present paper we do not present any 
parametrisation of the energy dependence of the value of the parameter \s{E), since 
amount of experimental data is not enough yet to check up this behaviour. 
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Figure 12: The experimental data on the rapidity dependence of the production denstities of A 
hyperon produced in proton collisions on Ar and Xe nuclei |73) (upper panel) and of A and A hyperons 
produced in proton collisions on S and Au nuclei [TH [75l [76] (lower panel), all at 200 GeV/c, together 
with the results of the corresponding QGSM calculations (see the main text por the description of the 
different curves). 
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